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The first step in the building of a spacetime solution of Einstein's gravitational 
field equations via the initial value formulation is finding a solution of the Einstein 
constraint equations. We recall the conformal method for constructing solutions 
of the constraints and we recall what it tells us about the parameterization of 
the space of such solutions. One would like to know how to construct solutions 
which model particular physical phenomena. One useful step towards this goal 
is learning how to glue together known solutions of the constraint equations. We 
discuss recent results concerning such gluing. 



1 Introduction 

The constraint equations of Einstein's theory of the gravitational field are 
familiar to most of those who work in general relativity: If we choose a set of 
appropriate local spacetime coordinates (x a , t) and if we write Einstein's field 
equations for the gravitational metric field g together with a coupled matter 
source field^] $ in the tensorial formf] 



then the four equations 



G F (9) = /cV( 9i f), (1) 



G tv {g) = KT tv {g^), (2) 



involve no time derivatives of higher than first order. These are the Einstein 
constraint equations. 

What makes the constraint equations important in the study of the 
physics of the gravitational field is their role in the initial value formulation of 



a We presume that the matter source field is not derivative-coupled, and that the coupled 
field equations sanstitute a well-posed PDE system. 

b We use MTWlH conventions on signs of the curvature and on indices: Greek letters run 
from to 3 (spacetime) while Latin letters run from 1 to 3 (space). Here n is a coupling 
constant. 
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Einstein's equations. We recall how the initial value formulation works: Say 
we want to construct a spacetime solution {g,^} of the Einstein equations 
(|l|) on a manifold M 4 = E 3 x R, where E 3 is a smooth three-dimensional 
manifold. The first step in doing this is to find a set of initial data consisting 
of 

jab - a Riemannian metric on E 3 

K cd - a symmetric tensor field on E 3 

ip - space covariant pieces of the matter field \P on E 3 
with {7, K, t/j} satisfying the constraint equations (0). In terms of the data 
{7, K, -0} the constraint equations take the form 

R + {trKf - K cd K cd = 2^(7, 0) (3) 

VaKS - V b {trK) = k J 6 ( 7 , 0) (4) 

where V a is the Levi-Civita^ covariant derivative corresponding to 7, R is its 
scalar curvature, p is a scalar function of 7 and if>, and Jb is a vector field 
function of 7 and -0-0 

Once we have obtained constraint-satisfying initial data {7, K, ip}, we may 
evolve this data in time by choosing everywhere on E 3 a scalar field N (the 
"lapse function" ) and a vector field M a (the "shift vector" ) and then solving 
the Einstein evolution equations in the form 

^■7q6 = ~2NK ab + £ M 7ab (5) 

and 

±K% = N{R c d + trKKj) - V c V d N + nT c d + £ M K c d (6) 

where T d are the spatial components of the stress-energy tensor corresponding 
to the matter field We note that the time evolution of the lapse and shift 
are not determined by the Einstein equations; they may be chosen freely, and 
they must be chosen for all t if one wishes to evolve 7 and K in time. The role 
of N and M is to determine the coordinates relative to the evolving spacetime. 

c i.e., metric compatible and torsion- free 

d \{, for example, the "matter source field" is Maxwell's electromagnetic field, then the 
spacetime covariant field is the spacetime one- form field A u , the space covariant fields 
1/1 may be chosen to be the magnetic spatial one-form field B a and the electric spatial 
vector field E a , and one calculates p = \(E a E b 'y ab + B a B b -y ab ) and J c = r] cmn E m B n , 
where rj cmn is the alternating symbol 77123 = 1, etc. In addition, the data must satisfy the 
Maxwell constraint equations SJ a B a = and V a E a = 0. 
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After N(x,t) and M c (x, t) have been chosen and the evolved data jabix, t) 
and K cd {x,t) and ip(x,t) have been determined^], we may reconstruct the 
spacetime metric (on the manifold M 4 = S 3 x /, where I is the interval for 
which the evolution can be carried out) via the formula 

g(x,t) = -N 2 dt 2 + j ab {dx a + M a dt)(dx b + M b dt). (7) 

We similarly reconstruct It then follows (from the Gauss-Codazzi-Mainardi 
3+1 decomposition equations for the spacetime curvature) that g and '5 satisfy 
the Einstein equation ([!]) on M 4 . 

Does the initial value formulation work? This general question is best 
considered by breaking it up into a number of more specific questions: 

1. Does every choice of initial data {j,K, -0} which satisfies the constraint 
equations evolve into a spacetime solution {g, ^}? 

2. How do we find data {7, K, tp} which satisfy the constraint equations 
©-(1? 

3. What is the space of solutions of the constraint equations? 

4. If we wish to use Einstein's equations to model the physical behavior of 
a particular physical situation (e.g., a black hole collision), how do we 
build that physical situation into the choice of initial data? 

5. What criteria should we use in making a choice of the lapse and shift? 

6. What do we know about the long time behavior of solutions? 

7. Is the initial value formulation a practical way to construct and study 
solutions of Einstein's equations? 

In this essay, we shall focus our attention on questions 2, 3, and 4, which 
concern finding and understanding solutions to the constraint equations. We 
start in Section 2 by discussing the conformal method for constructing so- 
lutions to the constraints. After recalling the basic steps of the conformal 
method, we state the results which tell us the extent to which the method is 
known to work. Based on these results, we discuss (also in Section 2) what is 
known about the function spaces which parameterize the set of solutions. 

Knowing how to construct constraint-satisfying initial data sets from con- 
formal data sets, and knowing how to parameterize the space of all such initial 
data sets, does not tell us much about how to use solutions constructed via the 

e The coupled Einstein-matter source field equations determine evolution equations for if). 
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initial value formulation for studying physical questions. With observational 
data from LIGO and other gravitational radiation detectors expected to start 
arriving within the next few years, it is crucial that we learn how to do this. 
Although progress to date on this front has been fairly limited, recent results 
on the gluing of initial data sets could be useful. We discuss some of these 
gluing results and the ideas involved in Section 3. 

Only if initial data can be evolved is it worthwhile to construct it. We 
comment briefly on evolution in Section 4, touching on questions 1, 5, 6, and 
7. We make some concluding remarks in Section 5. 



2 The Conformal Method and the Space of Solutions of the 
Constraint Equations 

For the past thirty years, the conformal method has been the most success- 
ful procedure for producing and studying solutions of the Einstein constraint 
equations, both numerically and analytically. The key idea of this method. 
which was developed primarily by Lichnerowicz, Choquet-Bruhat and Yorko 
is to split the initial data {7, K, ip} into two parts: The first part is to be cho- 
sen freely, while the second part is to be determined by the constraints. As 
a consequence of the split, the constraint equations take the form of a deter- 
mined PDE system of four equations for four unknowns (equations (|l2|)-(|l3]) 
below) , rather than the form of an underdetermined system of four equations 
for twelve unknowns (equations (||)-(f|)). 

To illustrate how this method works explicitly with matter source fields 
as well as gravitational fields, let us consider the Einstein-Maxwell equations, 
for which the initial data consists of {j a b, K cd , B a ,E a } and for which the full 
set of constraint equations takes the form (presuming no charges present) 

R + (trK) 2 - K cd K cd = n{B a B a + E a E a ) (8) 

V„JC? - V b {trK) = KVbmn B m E n (9) 

V a B a = (10) 

V a E a =0 (11) 
For this system, the split of the data is as follows: 
Free ("Conformal") Data: 

X a b - a Riemannian metric (specified up to conformal factor) 
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a cd - a divergence-free^] and trace-free symmetric tensor field 
(V c a cd = 0) (X cd ca cd = 0) 

t - a scalar field 

j3 a - a divergence-free vector field 
e a - a divergence-free vector field 

Determined Data 

4> - a positive definite scalar field 
T4^ a - a vector field 

Determining Equations 

V a (LW) a b = ^ 6 V b r + K Vbmn f] m e n (12) 

V 2 = - i(<7 a6 + LW^)(£T„6 + W afc )0~ 7 

-|(/? Q /3a+£ a ea)0- 3 + ^r 2 5 (13) 

Here the covariant derivative, the Laplacian, the scalar curvature, and the 
index manipulations all correspond to the metric X ab . Note that L is the 
conformal Killing operator, defined by 

LW ab := V a W b + V b W a - ^X ab y c W c . (14) 



The determining equations (|12[)-(|13[) are to be solved for cf> and W a . If, 
for a given set of conformal data {X ab ,cr cd ,t, (3 a 7 e a } one can indeed solve 
(|l2|)-(|l3|), then the reconstructed data 

lab = 4> 4 Xab (15) 

R cd _ + Lw cd^ + _^-A x cd T ^ 

3 

B a = 4>- 2 X ab p b (17) 

E a = <f>- 6 e a (18) 



■^In the conformal data, the divergence-free condition is defined using the Levi-Civita co- 
variant derivative compatible with the conformal metric A ct j. 
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constitute a solution of the constraint equations (|8|)-(|ll|). 

Does the conformal method always work in the sense that, for any choice of 
conformal data {X a b, o~ cd , t, f3 a , e a } one can always solve @-(|l3|) for {(j>, W a } 



and thereby obtain a solution {X a b, K cd , B a , E a } for the constraints (p|)-(pl|)? 
It is easy to see that this is not the case, as illustrated by the following simple 
example: 

Let us choose the manifold E 3 to be the three-dimensional sphere S 3 , 
and let us choose the conformal data to consist of the round sphere metric 
A = gr ro und) w ith R = I, together with <r cd = 0, e a = 0, f3 a = 0, and r = 1. 
With t constant and with r]bmn£ m P n equal to zero, the right hand side of 
equation ( |l2| ) vanishes. Since the operator V • L is self-adjoint and injective 
(up to conformal Killing vector fields) , it readily follows that the solutions W a 
to ([l2|) for this choice of conformal data all satisfy the condition LW a b = 
(which is the conformal Killing equation). Using this result, we find that the 
remaining equation ( |l3| ) takes the form 

V ^=^+l2^- ( 19 ) 

Since we seek solutions <f> which are positive definite, for any such solution 
the right hand side of equation ( |l9[ ) is positive define. But the maximum 
principle tells us that there are no functions on a compact manifold for 
which V 2 > (or < 0). Hence, for this choice of conformal data, there is 
no solution {(f>, W a } to equations (|l2|)-(|l3|), and therefore no corresponding 
solution to the Einstein-Maxwell constraint equations. 

This example shows that there are choices of conformal data which are 
not mapped to solutions of the constraint equations by the conformal method. 
What happens generally? This question is difficult to address in any complete 
sense because there is a very wide variety of cases of physical and mathematical 
interest, and each must be handled separately. The various cases may be 
classified using the following criteria: 

Manifold and Asymptotic Geometry 

• Closed E 3 

• Asymptotically Euclidean 

• Asymptotically hyperbolic 

• Compact E 3 with boundary conditions 

Mean Curvature 

• Maximal 
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• Constant non maximal (CMC) 

• Non Constant (non CMC) 



Fields 

• Einstein vacuum 

• Einstein - standard matter source field£] 
Function Spaces 

• C a (analytic) 

• C°° (smooth) 

• C k+fi (Holder) 

• Hi (Sobolev) 

• H% g (weighted Sobolev) 

For a number of cases (labeled using these criteria), we have a fairly 
complete understanding of which sets of conformal data map to solutions and 
and which do not. This is true for the following: 

1. Closed E 3 ; CMC; Einstein vacuum or Einstein-standard matter source 
field; A G C 3 and a, s, (3 G H$(p > 3). 

2. Asymptotically Euclidean; Maximal; Einstein vacuum or Einstein- 
standard matter source field; Data in weighted Sobolev spaces. 

3. Asymptotically hyperbolic: CMC; Einstein vacuum or Einstein-standard 
matter source field; Data in weighted Sobolev spaces. 

For these three cases, one finds that all sets of conformal data which are 
not readily disqualified by simple sign criteria or by the maximal principle are 
mapped to solutions by the conformal method. To explain the sense in which 
this holds, let us focus for the moment on the first of these cases; that with 
E 3 closed: 

The Yamabe classification of metric conformal classes is a key tool for de- 
termining which sets of conformal data on a closed manifold E 3 are mapped 

9 We label as "standard matter source field" any field theory which couples to Einstein's 
theory without derivative coupling. Included are Einstein-Maxwell, Einstein- Yang-Mills, 
Einstein-Dirac, Einstein-Klein-Gordon, Einstein-Cartan, and Einstein fluids. For all these 
theories, there are conformal splittingpi such that, if Vt = 0, the constraint equations 
semi-decouple. See Isenberg and NesteiQ. 
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to solutions of the constraints. This classification is based on the Yamabe 
TheoremH which says that every metric X a b on a closed manifold E 3 is con- 
formally related to another which has constant scalar curvature either +1, 
0, or —1; and moreover, the sign of R for constant scalar curvature metrics 
conformally related to a given X a b is unique. Thus the set of all metrics (and 
metric conformal classes) on a closed £ 3 is partitioned into the three Yamabe 
classes y + (Y, 3 ), 3^°(S 3 ), and y~(Y, 3 ) depending on this sign. 

The Yamabe theorem and Yamabe classification are important for 
CMC conformal data on a closed S 3 because solubility of (|l2|)-(fi"3|) 
is invariant under the conformal transformation {X a b, o~ cd , t, j3°, e a } — > 
{d i X ab ,e- lo a cd ,T,0- 6 f3 a ,e- 6 e a } for any positive 9. (This is not true if 
Vr ^ 0). Thus, to determine for which sets of conformal data (|l2|)-(|l3}) admit 
a solution, one may without loss of generality restrict attention to conformal 
metrics with R(X) — —1,0, or +1. 

As a consequence of this conformal invariance and the resulting sim- 
plification, one immediately sees from ( |l3| ) and the maximum principle 
that the following categories of conformal data do not lead to solutions^ 
{y + ;(a,p) = 0;r = 0}, {y-\(a,p) = 0;r = 0}, {y+;(a,p) = 0; r ? 0}, 
{y°;(<j,p) =0;r^0}, {y°; (a, p) £ 0; r = 0}, {y~; (a, p) # 0; r = 0}. Much 
less evident is the fact that for all other sets of conformal data on closed S 3 
(|l2|)-(|l3|) can be solved, and we obtain a solution of the constraints. The 
proof of this result (with one. case left out) is discussed in the review paper 
of Chofluet-Bruhat and YorkB; the complete proof appears in the work of the 
authorB. 

Analogous results, with a few important modifications, hold for the other 
two cases mentioned above. Specifically, for the asymptotically Euclidean case 
with t = 0, we first note that there is a Yamabe-type partition of asymptoti- 
cally Euclidean metrics into two classes Z + and Z~ , with the Z + conformal 
geometries admitting a conformal transformation to metrics with R = 0, and 
the Z~ conformal geometries nx>t admitting such a transformation. It follows 
from work of Brill and CantorQ that asymptotically Euclidean metrics with 
non-negative scalar curvature, as well as those satisfying a cextain integral 
condition, are contained in Z + . The work of CantorQ (See alsota) shows that 
a set of asymptotically Euclidean conformal data {X a b, °~ cd , t, (3 a , e a } (in ap- 
propriate function spaces) is mapped to a solution by the conformal method 
if and only if A G Z + (The behavior of the other fields is irrelevant.). 

For the asymptotically hyperbolic case (with t ^ 0), there are important 



h Here (<r, p) = means that a cd and p are identically zero everywhere on S 3 , while (tr, p) =a= 
means that there exists at least one point on S 3 at which cr cd or p is not zero. 
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questions regarding the asymptotic behavior of solutions; but ignoring those, 
the results are very simple: For every set of asymptotically hyperbolic confor- 
mal data, emulations ([E|)-(|l3|) can be solved, and we obtain a solution to the 
constraintsB LL9. 

These results together tell us that if we are interested in constant mean 
curvature (or maximal) initial data, the conformal method works, apart from 
a few readily identified special cases. This is useful for two reasons. The first 
is that since equations of the form of the determining equations (p^|)-(p^|) can 
be handled by modern computer algorithms (at least when they are semi- 
decoupled, which holds so long as we restrict to conformal data with constant 
t), the conformal method is an effective way to produce CMC initial data sets 
for numerical relativity. The second is that, as a consequence of uniqueness 
theorems relating conformal data and solutions of the constraintsO, together 
with the results noted above, certain function spaces of conformal data pro- 
vide an effective parameterization of the set of solutions of the constraintsEil. 
In light of the uniqueness of CMC foliations of spacetime solutions of Ein- 
stein's equationsEj, one may use these function spaces of conformal data to 
parameterize the set of spacetime solutions as welL/ignoring those solutions 

10). 



which do not admit CMC or maximal foliationsEj E_l) . 

These results are very useful, but they are limited in a number of ways, 
and so there is a lot more that needs to be done. First, one would like to 
know which non CMC conformal data sets lead to solutions. While there are 
some results concerning this question, they all presume that |Vr| is small. 
Indeed, for small |Vr|(with certain restrictions on the zeros of t1 
show that solutions to the coupled system (|l^)-(13) always existt 
Unfortunately, the iteration methods we have used to prove the results just 
cited do not appear to be able to handle general |Vr|, so new ideas are likely 
needed. 

Next, we would like to be able to generate solutions to the constraint 
equations on manifolds with boundaries. The initial value-boundary value 
problem has been studied for Einstein's theory in Friedrich's formulation^, 
but there has not been much mathematical work on it in the standard {7, K} 
formulation discussed here. In view of the recent numerical work on space- 
times with excised black holes, a mathematical treatment, in terms of the 
{7,-ftT} formulation, of solutions with boundaries would be useful. 

A third issue regarding the construction of solutions of the constraint 
equations which should be addressed further is the loosening of differentiabil- 
ity requirements for solutions, including weak solutions. The motivation for 
doing this comes_fxom the progress that has recently been made by Klainerman 
and Rodnianskio and by Smith and TataruEJ in proving local well-posedness 



iff 
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for the Einstein equations for data of lower differentiability. The goal (based 
on energy considerations) is H^, 2 ; so far -fff+e nas been obtained. This leads 
one to tryto construct initial data sets in iff+e or even H^ 2 - The results 
to dateS i produce initial data sets which are considerably more restricted 
(i.e., smoother) than this, but work is proceeding toward obtaining solutions 
of the constraint equations with this low differentiability. 

The question that remains is how to construct initial data sets which can 
be used to model physical phenomena of interest. We discuss this in the next 
section. 



3 Physical Modeling and the Gluing of Solutions of the 
Constraint Equations 

Say we want to study the gravitational radiation produced by the inspiral 
collision of a neutron star and a black hole. The first step of a numerical 
modeling of this phenomenon is to find initial data {7, K, ip} which constitutes 
a snapshot of the two objects and their ambient spacetime at some moment 
prior to the collision. It is important that the initial data accurately represent 
the two objects and their surroundings in a relatively quiet pre-collision state, 
without any extraneous unrelated radiation, or the modeling will be of little 
use. 

The electromagnetic analogue of this modeling problem is familiar and 
easily handled, since the background space is fixed and independent of the 
electromagnetic fields in Newtonian-Maxwell theory, and since the represen- 
tation of electromagnetic radiation is well understood in this theory as well. 
So, too, if we work with a post-Newtonian approximation to Einstein's theory, 
the choice of ini±iaLdata for phenomena like the inspiral collision is understood 
to some extentE3 Ei 

However, attempts to find physically accurate initial data sets for phe- 
nomenon such as this in terms of the full Einstein theory have been stymied, 
despite many years of significant effort, both numerical and analytical. The 
conformal method can of course be used to generate candidate sets of data, 
and such sets are used in most numerical modeling studies. The lack of a priori 
control of the conformal factor and therefore of the ambient space, together 
with the difficulty in representing and recognizing gravitational radiation in 
the gravitational radiation in the conformal data {A, a, t, if)} for massive rel- 
ativistic objects renders these sets somewhat suspect, however. 

One approach that has been discussed to try to control the choice of 
initial data is to work with sets of post-Newtonian initial data in the very 
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early pre-radiation stages of the collision and then try to evolve these into the 
later, more relativistic stage. This approach has some promise, but can be 
very expensive in numerical evolution time. 

Another possible approach that could be tried relies on "gluing." The idea 
here is to start with well-understood initial data sets for separate portions of 
the spacetime-say, one for the black hole and one for the neutron star-and 
develop a procedure for joining these together into one set. Mathematically, 
the idea is to join two or more solutions of the constraint equations into a 
single solution in such a way that, away from the joining region, the data is 
largely unchanged. Physically, the challenge is to understand and control the 
physical effects of the joining region. 

There has been substantial progress during the past couple of years on the 
mathematical aspects of gluing solutions of the Einstein constraint equations. 
Mazzeo, Pollack, and the author (IMP) have proven theorema_which prescribe 
a procedureJbt-gluing connected sums of sets of initial dataE3, while Corvino 
and SchoerO E3 have shown that rather general asymptotically flat initial 
data sets can be glued to exact Schwarzschild or exact Kerr exteriors outside 
a transition region. We first discuss the IMP work, and then briefly describe 
the results of Corvino and Schoen. 

3.1 Connected Sum Gluing 

We start with a pair of constant mean curvature solutions {E 3 , 71, ifi} and 
{£2, 72, K2} of the vacuum constraint equations (3)-(4). (We presume for now 
that there are no matter source fields tp, and so p — and J = 0). Each of 
these solutions may be asymptotically Euclidean, asymptotically hyperbolic, 
or one may have either or both of £ 3 and Ef, closed. The two solutions need 
not have matching asymptotic properties, but they do need to have the same 
constant mean curvature r. 

If we now pick a pair of points pi € £ 3 , P2 € £2 on each manifold, then the 
gluing procedure produces a one-parameter family of solutions {E 3 ,7 S , K s } 
with the following properties: (i) E 3 = Ej^E 3 ,, with £ 3 joined to £2 by a 
"neck" (S 2 x R) connecting a neighborhood of p\ to a neighborhood of pi 
(both now excised); (ii) the data K s } can be made to be arbitrarily close 
to the original data {71, K\] and {72, K2} on appropriate regions of £ 3 away 
from the neck by choosing s sufficiently large; (iii) trK = r; and (iv) the 
geometry on the neck (i.e., lengths, curvature) can degenerate for s — > 00, but 
its behavior is controlled by exponential bounds in s. 

This gluing procedure works for quite general classes of data. The only 
restrictions on {£ 3 , 71,^1} and {£3,72,^2} besides the CMC matching con- 
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dition trK% = trK2 and. certain regularity restrictions (stated in terms of 
weighted Holder spaces Ea) are that if either of the sets of data has S 3 closed, 
that set of data must have K cd not vanishing identically, and must be free of 
nontrivial conformal Killing vector fields|. Neither of these latter conditions 
are needed for asymptotically Euclidean or asymptotically hyperbolic data, 
and we note that if a set of data {£ 3 ,7, K} fails to satisfy either condition, 
it will generally satisfy them both after a small perturbation is made. 

The procedure allows a number of interesting sets of initial data to be 
produced. One can, for example, add a sequence of small black holes to 
any asymptotically flat spacetime {M 4 ,<?} by gluing a sequence of copies 
of Minkowski data to asymptotically Euclidean data for {M 4 ,g} and then 
evolving the glued data. These_spacetimes (somewhat reminiscent of the 
Misner multi-black hole solutionsE2l, but far more general) may be constructed 
with the interiors either independent or connected. Another set of spacetimes 
one obtains via this gluing procedure allows one to study the question of black 
holes in cosmological spacetimes: One can glue asymptotically Euclidean or 
asymptotically hyperbolic data to cosmological data on closed £ 3 and then 
see if black hole-like physics develops^]. This gluing procedure also permits 
us to add one or more wormholes to most solutions of Einstein's equations: 
If one chooses the pair of points p\ and P2 on a single connected initial data 
set {S 3 ,7,iC} and carries out the gluing procedure, one obtains initial data 
much like {T, 3 ,j,K} except with an S 2 x R wormhole glued on. 

To what extent can one carry out these gluings in practice, say, numeri- 
cally? To answer this question, we shall now describe the basic steps of the 
connected sum gluing procedure: 

Step 1: One starts by performing a conformal transformation of the data 
lab —* "0~ 4 7a6, o cd — ► ijj 10 a cd , (Here a cd is the trace-free part of K cd .) which 
is trivial (ip — 1) away from a neighborhood of each of the points pi and p2, 
and which is singular at each of the points (with a specified type of singular 
behavior). The effect of this conformal blowup is to replace a ball surrounding 
Pi with an S 2 x R infinite length half-cylinder which is asymptotically a 
standard round half- cylinder, and to do the same in a neighborhood of pi- 
Note that there are standard explicit formulas one can choose for ip to carry 
out this first step. 

Step 2: Each of the newly added half cylinders extends out infinitely 
far. If one cuts each at a distance s/2 (measured using the local conformally 

! Actually, a Killing vector field can be present, so long as it does not vanish at the chosen 
points pi or p2- 

i Note that there is no accepted definition of a black hole in a spacetime without some sort 
"asymptotic infinity" structure. 
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transformed metric) out from the neighborhoods of pi and P2, and identifies 
the two cut half cylinders on the joining sphere, then one has the manifold 
E 3 = E 3 ^!] 3 ,, with its connecting s-parametrized "neck". Then, using a pair 
of cut-off functions xi an d X2i one may patch together the metric -% = 71X1 + 
72X2 and the traceless apart of the second fundamental form & s = o-\X\+o~iXi- 
Both 7 S and a s arc well-defined on E 3 , and both depend on the gluing distance 
parameter s. One also glues together the conformal function ^ S p. i/'iCi+V^^, 
using a different set of cut-off functions. (For details, see IMPc5). This step, 
like the first, can be done using standard explicit formulas for XiiX2,£i, and 

6- 

Step 3: A careful choice of the cut-off functions results in a s be- 
ing everywhere traceless with respect to The conformal transformation 
<j cd — > -0 lo cr cd guarantees that away from the neck, <x s is divergence-free with 
respect to However, in the neck, V c af 7^ 0. So one replaces of by (if, 
which is obtained by solving 

V c (LY s ) cd = V c of (20) 

for the vector field Yf, and setting 

(xf = of - (LY s ) cd . (21) 

The tensor field (if is divergence-free everywhere on E 3 with respect to j s . 
This step is straightforward to carry out in practice. However, to be able to 
proceed further with the gluing procedure, one needs to know not only that 
the divergence-free field (if exists, but also that || (i s — a s \\c k +e is very small 
for sufficiently large s. It is shown in IMP that this estimate always holds. 

Step 4: {j s ,(i s ,t} is a standard set of conformal data on S 3 . So one 
may set up and attempt to solve 

V 2 s = ±R.4> a - ^(i scd (if<P~ 7 + ^r 2 ^ (22) 

for cj) s . Although the Yamabe class of the metric % constructed in Step 2 is 
not evident, the estimates obtained in Step 3 allow one to show that (|2^) has 
a unique solution. Further, with a considerable amount of work, one shows 
that || <j) a — ip s \\c k +P is very small for sufficiently large s. As noted above, 
once one knows that a solution to ( p2| ) exists, one can readily obtain that 
solution numerically. 

Step 5: Using formulas analogous to (p^|)-(^6|), one constructs 7, and 
K s from <j) s and {7 s ,/t s ,r}. One immediately verifies that for each value of 
the parameter s, {E 3 , 7 S , -ftT s } solves the vacuum constraint equations. As a 
consequence of the estimates of Steps 3 and 4, one verifies that away from the 
neck, 7, and K a approach the original data for sufficently large s. 
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We see from the description of these five steps that in practice, this con- 
nected sum gluing procedure involves making a few choices of standard cut 
off functions, plus solving familiar elliptic equations. It should be a useful 
procedure for studying solutions of Einstein's equations numerically. 

We note that, to date, the IMP connected sum gluing procedure has been 
proven to work for the vacuum Einstein case only. It is expected that it 
can be implemented for Einstein-Maxwell, Einstein- Yang- Mills, and the other 
Einstein-standard matter source fields as well; work is under way to show this. 
It would be nice to also extend the procedure to non CMC data, but this is 
likely to be difficult. 

3.2 Exterior Schwarzs child and Kerr Gluing 

The starting point for the exterior Schwarzschild gluing procedure is any 
asymptotically Euclidean, time symmetric (K cd = 0) initial-data set {E 3 ,7}. 
One chooses a compact region A 3 in E 3 . One then shows E£l that there is an 
asymptotically Euclidean time symmetric solution {E 3 , 7} such that (i) inside 
a ball B r which contains A 3 as a subset, 7^ = 7; and (ii) outside the ball 
B2r, we have 7|s 3 \B 2 r = lsch{m), where ~/s c h(m) is the spatial Schwarzschild 
metric for some positive mass m. 

We note three important features of this type of gluing: First, when the 
gluing is completed, the region inside B r as well as the exterior region outside 
Bir are unchanged. This is not the case for the connected sum gluing, in which 
the original solutions are slightly changed (away from the neck) by the gluing 
procedure. Second, the analysis used to prove this result does not rely in any 
way on the conformal method. That is, one does not construct the solution 
{E 3 ,7} by first specifying a conformal metric everywhere and then solving 
for the conformal factor. Indeed, one can carry out exterior Schwarzschild 
gluing only if one works with the constraints as an underdetermined system 
to be solved for 7, rather than as a determined system to be solved for 4>. 
Third, the proof that this gluing works does not readily translate into a step 
by step procedure that one can carry out numerically. The proof guarantees 
existence, but does not prescribe a procedure for constructing {E 3 ,7} 

The result is quite surprising and remarkable. It says that for any com- 
pact portion of a time symmetric asymptotically Euclidean solution, one can 
arrange the gravitational fields in an annular region around it so that outside 
this region, the gravitational field is exactly Schwarzschild. Note that the 
analogous result does not hold true for Newtonian gravity. Note also that, 
combining this result with the work of Friedrichcj, one may be able to pro- 
duce a very large family of solutions of Einstein's equations with a complete 



GR16Final: submitted to World Scientific on February 7, 2008 14 



X (null asymptotic infinity). 

One large but unsurprising restriction on the type of data which can be 
glued to an exterior Schwarzschild spatial slice is that the data have K cd — 
0. More recent results of Corvino and Schoen allow one to remove this 
restriction provided that one replaces the Schwarzschild exterior by a slice of 
the more general Kerr (rotating black hole) They do impose restrictions on 
the data {£ 3 ,7,AT} which can be patched to an exterior Kerr slice: They 
require certain decay behavior in 7 and K which are essentially just enough 
so that the momentum and angular momentum integrals at infinity are well- 
defined. These restrictions appear to be minor, in view of the result. Again, 
one consequence of these results might be an even larger family of solutions 
with complete I. 

4 Some Comments on Evolution 

Of the seven questions raised in the introduction, we have focussed in this 
essay on those three - 2,3, and 4 - which deal with the constraint equations. 
We shall now briefly comment on questions 1,5, and 6, which mostly concern 
evolving sets of initial data into spacetime solutions. 

The first question addresses the issue of local existence and well-posedness 
of the initial value formulation of the Einstein equations^ To a large extent, 
this issue was settled fifty years ago by Choquet-BruhaiE3 who proved (using 
harmonic coordinates) that for smooth initial data (C 1 , with large I) the 
Einstein system is well-posed. While this result is sufficient for most purposes, 
there have been strong attempts in recent years to prove well-posedness for 
more general classes of data. The goal is to prove it for data-in the Sobolev 
space H \ . The recent work of Klainerman and RodnianskiE3, and of Smith 

and Tataru0achieves it for -ff| +£ (e > 0) data. 

Well-posedness theorems show (among other things) that data in the spec- 
ified function space can always be evolved to a spacetime solution for a suffi- 
ciently small time duration into the future and into the past. They say little 
about how long a spacetime solution will last, or what its behavior will be 
far into the future and far into the past. While these issues are very difficult 
and far from resolved, there has been substantial progress on them in recent 
years. We note in particular the work of Christodoulou and KlainermanEll, 
which verifies the stability of Minkowski spacetime (and thereby proves that 
there is a nontrivial family of asymptotically flat spacetimes which extend 
for infinite proper_time into the future and into the past); that of Anders- 
son and Moncrie£3 which does the same for the stability of the expanding 
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k = — 1 vacuum Friedman.- Robertson- Walker cosmological spacetimes; and 
the work of ChristodoulouEa which verifies weak cosmic censorship^ for spheri- 
cal collapsing Einstein-scalar field spacetimes. There has also been a collection 
of works by Berger, Chrusciel, Garfinklc, Kichcnassamy, Moncrief, Rendall, 
Ringstrom, Wainwright, Weaver, and the author which provide increasing ev- 
idence for the presence of asymptotically velocity dominated behavior and 
mixmaster behavior in cosmological spacetimes, and further tends to support 
the validity of the strong cosmic censorship coajectureH in these spacetimes. 



The construction or study of spacetime solutions of Einstein's equations 
via the initial value formulation requires one to make a choice of spacetime 
foliation (i.e., a specification of the t — constant hypersurfaces which fill the 
spacetime) and of spacetime threading (i.e., a specification of the x a — con- 
stant observer world lines which are everywhere transverse to the foliation, 
and fill the spacetime). The choice of foliation is controlled infinitesmally by 
the lapse function N, and the choice of threading is similarly controlled by 
the shift vector M a . 

In constructing a spacetime from specified initial data, one may always 
make the simple choice N = 1 and M a — ("Gaussian normal coordinates"). 
However, this choice generally leads to a non physical and premature break- 
down in the pplution (a coordinate singularity) . Consequently, the maximal 
developmentEZl of the initial data may not be obtained. So, one criterion for a 
good choice of the lapse and shift is that the choice effectively avoids such co- 
ordinate singularities. Also important is that the choice be relatively easy to 
implement in practice, and that it not obscure the gravitational physics of the 
spacetime (e.g., by simulating gravitational radiation which is not physically 
present). 

Constant mean curvature or maximal slicing, coupled with some sort of 
coordinate shear minimizer, is often cited as a choice of foliation and thread- 
ing which avoids coordinate singularities and clarifies the physics. However, 
implementing this choice requires that one solve a set of elliptic equations 
on each time slice. This is a large expense in computer time for numeri- 
cal constructions, and it precludes explicit forms for the lapse and shift in 
analytical studies. For special families of solutions such as the Gowdy space- 



fc The weak cosmic censorship conjecture says that in asymptotically Euclidean spacetime so- 
lutions, the singularities which develop during gravitational collapse are generically shielded 
from the view of observers at infinity by tkc development of an event horizon. 
'The strong cosmic censorship coniecturePH says that in generic solutions developed from 
Cauchy data on a compact Cauchy surface, a Cauchy horizon (with its attendant causality 
difficulties) does not form. 
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times or the U(l) Symmetric solutions , there are certain choices of lapse and 
shift picked-out by the geometry (areal for Gowdyll, harmonic time for U(l) 
SymmetricE3) ; but these are special cases. More generally, the choice of lapse 
and shift remains a difficult issue. 

5 Concluding Remarks 

We have discussed a number of the challenges that one encounters in con- 
structing and studying spacetimes via the initial value formulation. These 
occur both in finding sets of initial data which satisfy the Einstein constraint 
equation, and in evolving these sets of data. Some are fairly difficult. How- 
ever, paraphrasing Winston Churchill's description of democracy, one finds 
that the initial value formulation is the most impractical way to work with 
solutions of Einstein's equations, except all of those other ways which have 
been tried from time to timeo. The fact that most numerical studies of solu- 
tions are carried out using the initial value formulation attests to the relative 
practicality of this approach. 
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